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ABSTRACT 

We  present  a  convenient  method  for  determining  the  thickness  of  a  cryogenic 
parahydrogen  (pH,)  solid  from  its  infrared  (IR)  absorption  spectrum.  Millimeters-thick  pH, 
solids  of  exceptional  optical  clarity  can  be  produced  by  the  rapid  vapor  deposition  method  [M.E. 
Fajardo  and  S.  Tam,  J.  Chem.  Phys.  108,  4237  (1998)].  Doping  of  these  pH,  solids  is  readily 
accomplished  by  co-deposition  of  the  desired  impurities,  making  them  excellent  hosts  for  high- 
resolution  matrix  isolation  spectroscopy.  The  intensities  of  the  IR  “double”  transitions 
Q,(O)+S0(O)  and  S,(0)+S0(0)  of  the  pH,  host  are  insensitive  to  the  matrix  microstructure  and  to 
the  presence  of  dopants,  so  these  absorptions  are  especially  well  suited  for  thickness 
determinations.  We  calibrate  the  integrated  absorption  intensities  of  these  two  bands  against  the 
sample  thicknesses  determined  from  transmission  interference  fringes  appearing  in  the  same 
experimental  spectra;  we  report:  a[Qi(O)+S0(O)]  =  4.84(±0.13)xl0  14  cm7s,  and  a[S](0)+S0(0)] 
=  0.35(±0.02)xl0'14  cm3/s  (95%  confidence).  We  also  discuss  several  other  advantages  of  rapid 
vapor  deposited  pH,  solids  as  hosts  for  quantitative  IR  absorption  spectroscopy  of  dopant 


species. 
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I.  INTRODUCTION 

We  recently  reported  the  production  of  millimeters-thick  parahydrogen  (pH2)  solids  of 
exceptional  optical  clarity  by  rapid  vapor  deposition  of  pre-cooled  pH2  gas  onto  a  liquid  helium 
(lHe)  cooled  substrate.1,2  Our  sample-in-vacuum  geometry  permits  the  facile  doping  of  these  pH2 
solids  with  a  variety  of  atomic  and  molecular,  charged  and  neutral  species1  which  we  interrogate 
via  matrix  isolation  spectroscopy  (MIS).  This  work  is  part  of  the  U.S.  Air  Force's  High  Energy 
Density  Matter  (HEDM)  program;9  our  goal  is  to  demonstrate  practical  chemical  energy  storage 
in  cryogenic  solids,  and  ultimately  to  use  these  materials  as  advanced  rocket  propellants. 

Of  more  general  interest  are  the  advantages  of  using  solid  pH2  as  a  host  for  MIS  over  the 
traditionally  employed  rare  gas  solids  (RGS).  High  resolution  (Av/v  <  10 6)  infrared  (IR) 
vibrational  absorption  spectroscopy  in  solid  pH,  was  first  demonstrated  by  Oka  and  co-workers 
at  the  U.  of  Chicago  in  the  1980's,10  and  further  developed  in  collaboration  with  the  group  of 
Shida  and  Momose  at  Kyoto  U.11,12  The  relatively  "soft"  nature  of  the  pH2  host10,12  and  the 
typically  weaker  guest-host  interactions,  reduce  the  occurrence  and  importance  of  trapping  site 
inhomogeneities,  resulting  in  the  observed  sharp  absorptions.  The  correspondingly  reduced  cage 
effect  in  solid  pH2  allows  for  more  efficient  production  of  novel  chemical  species  by  in-situ 
photolysis  or  radiolysis, 12,13  and  suggests  a  myriad  of  possibilities  for  studying  condensed  phase 
chemical  processes  in  a  minimally  perturbing  environment.  Furthermore,  foreign  dopant  species 
induce  weak  but  measurable  IR  activity  in  the  neighboring  pH2  molecules;  thus,  in  addition  to  the 
direct  IR  absorptions  of  the  dopant  "solute,"  a  single  IR  spectrum  contains  information  about  the 
response  of  the  pH,  "solvent,"  as  well.1  Growing  recognition  of  these  advantages  is  resulting  in 
an  increasing  popularity  of  solid  pH2  as  a  host  for  MIS  studies. 
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Another  advantage  of  utilizing  solid  pH2  as  a  host  for  MIS  that  has  received  little  explicit 
mention  to  date  is  its  special  suitability  for  quantitative  absorption  spectroscopy.  By  this  we 
mean  the  facile  application  of  Beer’s  Law14  to  determine  absolute  dopant  concentrations  from 
experimental  absorption  spectra.  This  is  of  obvious  importance  to  our  HEDM  chemical  energy 
storage  effort,  as  well  as  to  other  analytical  applications  of  MIS.15  For  experimental  decadic 
absorbance  spectra  recorded  as  a  function  of  wavenumber,  v  (cm  '): 


A(v)  =-log10 


I(v) 

I0(v) 


(1) 


we  state  Beer’s  Law  as: 


A(v)  =  2l03  Cxd  (2) 

in  which:  ax(v)  is  a  molecular  absorption  coefficient  (cmr/mol),  Cx  is  the  concentration 
(mol/cm3)  of  species  X,  and  d  is  the  sample  thickness  (cm).  Alternative  definitions  and  units 
systems  are  also  in  common  usage.16,17  For  spectrally  isolated  absorption  bands,  the  dopant 
concentrations  are  most  accurately  calculated  via  the  related  expression: 


Cv  = 


2.303  f  A  (v)  dv 

_ Jband _ 

d  La*(i7)  di7 


(3). 


We  also  routinely  report  fractional  dopant  concentrations  normalized  by: 


X%  =  100% 


V  ^  hosl  j 


(4) 


or 


X  (ppm)  =  106 


V  ^host  J 


(5) 


J 


i 


P 


in  which:  NA  =  6.022xl023  #/mol  is  Avogadro’s  constant,  and  Nh0SI  is  the  number  density  of  the 
matrix  host  (NpH2  =  2.60x1 022  #/cm3  at  lHe  temperatures,  ref.  1 8). 

Several  favorable  properties  of  rapid  vapor  deposited  pH2  solids  contribute  to  the 
straightforward  evaluation  and  accurate  application  of  these  expressions.  Their  extreme  optical 
clarity  eliminates  troublesome  sloping  spectral  baselines  usually  encountered  in  thick  rare  gas 
matrices  due  to  wavelength  dependent  scattering  losses.  The  (often)  sharp  dopant  spectral 
features  further  facilitate  the  determination  of  absorption  band  integration  limits  and  baselines. 
The  millimeters-long  useful  optical  pathlengths  and  sharp  dopant  spectral  features  combine  to 
permit  the  facile  detection  of  many  dopants  at  very  low  (~  1  ppm)  concentrations,  which 
eliminates  “deviations”  from  Beer’s  Law  due  to  dopant-dopant  interactions  (conversely,  dopant 
aggregation  and  other  “intermolecular  interactions”  can  be  studied  systematically  by  increasing 
concentrations  from  this  well  isolated  limit).  Analysis  of  the  condensed-phase  induced  IR 
absorptions  of  the  host  pH,  molecules19'28  demonstrates  that  the  as-deposited  pH2  hosts  are  close- 
packed  solids  with  a  mixed  hexagonal  close-packed  (hep)  and  face-centered  cubic  (fee) 
polycrystalline  microstructure,1,4  and  that  the  metastable  fee  regions  largely  anneal  away1,4'3,7,8 
upon  warming  above  T  ~  4.5  K.  Thus,  these  samples  are  neither  amorphous  nor  porous;  instead, 
they  have  very  nearly  the  same  bulk  density  as  pH2  solids  grown  from  the  melt.  As  we  will  show 
below,  the  pH2  absorptions  also  conveniently  permit  determination  of  the  sample  thickness  from 
the  same  spectral  data  containing  the  dopant  absorption  features,  ensuring  that  they  relate  to  the 
same  experimental  interrogation  volume.  Finally,  the  low  index  of  refraction  of  solid  pH2 
implies  only  weak  perturbations  to  the  absolute  intensities  of  dopant  vibrational  transitions,29'37 
reducing  the  magnitude  of  the  errors  associated  with  this  correction.  This  last  point  deserves  a 
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brief  farther  discussion. 

For  IR  absorptions  occurring  in  condensed  phases,  the  electromagnetic  field  experienced 
by  the  absorbing  molecule  is  modified  by  the  dielectric  properties  of  the  host.  Many  models  of 
the  consequences  of  this  “effective”  or  “local”  field  to  experimental  observables  have  been 
proposed.29'37  In  particular,  the  so-called  “Chako”  or  “Polo-Wilson”  equation29-30  can  be  used  to 
correct  the  molecular  absorption  coefficients: 


CCy 


a&  — 
«x 

n 


V+2^ 


(6) 


in  which  n  is  the  index  of  refraction  of  the  matrix  host.  Disagreement  persists  to  date  on  which 
“index  of  refraction”  should  be  used  in  Eqn.  (6);  em:  should  it  include  any  frequency 
dependence,  or  contributions  from  the  dopant  molecule,  etc.  For  our  purposes  we  take: 

n  =  no  =  ^  (7) 


in  which:  and  e0  are  the  DC  (zero  frequency)  refractive  index  and  dielectric  constant  of  the 

matrix  host,  respectively.  For  solid  pH,  at  lHe  temperatures  s0  =  1 .297  (refs.  38-40);  thus  = 

1 .14,  and  af2  =  1.06  .  For  comparison,  the  refractive  index  of  a  solid  Ar  matrix  is  n^  = 

1.29  (ref.  41),  resulting  in  a  significantly  larger  gas-to-matrix  correction  of:  =1.16  ajf* . 

In  what  follows  we  describe  our  approach  for  determining  the  thickness  of  a  pH2  solid 
from  the  condensed-phase  induced  IR  absorptions  of  the  host  pH,  molecules.  We  will  present 
data  gathered  from  MIS  experiments  on  over  100  different  pure  and  doped  pH2  samples.  We  will 
show  that  the  so-called  “double”  Q,(0)+S0(0)  and  S,(0)+S0(0)  transitions,42  in  which  a  pair  of 
nearest  neighbor  pH,  molecules  cooperate  to  absorb  a  single  photon,  are  well  suited  for 
thickness  determinations.  We  point  out  possible  pitfalls  to  using  “single”  transitions,  such  as  the 
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S,(0)  or  U0(O)  bands,  for  this  purpose.  We  provide  recommended  spectral  integration  regions  and 
baselines,  and  corresponding  calibration  factors  based  on  a  comparison  with  interferometncally 
determined  sample  thicknesses.  We  conclude  with  example  concentration  calculations  for  a  C02 

doped  pH2  solid. 

II.  EXPERIMENTAL 

Our  experimental  apparatus  and  techniques  have  been  described  before,  here  we  give 
only  a  brief  review  of  our  previous  descriptions,  and  focus  instead  on  aspects  relevant  to 
quantitative  spectroscopic  measurements. 

A.  Apparatus 

Fig.  1  shows  a  section  view  through  the  center  of  the  deposition  chamber  of  our  lHe  bath 
cryostat.  The  octagonal  vacuum  shroud  measures  7.62  cm  across  the  flats;  eight  o-ring  sealed 
ports  allow  for  sample  preparation  and  optical  access.  The  chamber  sits  directly  atop  a  small  (30 
1/s)  turbomolecular  pump  which  provides  a  room  temperature  vacuum  ~  10'7  Torr. 

Solid  pH2  samples  are  prepared  by  deposition  onto  a  0.4  cm-thick,  0.5°-wedged  BaF2 
substrate  clamped  in  a  gold-plated  oxygen  free  high  conductivity  (OFHC)  copper  holder  which  is 
cooled  by  the  lHe  reservoir;  all  thermal  connections  are  made  with  100  pm-thick  indium  foil. 

The  deposition  substrate  is  surrounded  by  a  liquid  nitrogen  (1N2)  cooled  aluminum  radiation 
shield  in  which  openings  are  cut  for  sample  preparation  and  optical  access.  Operating  the 
ortho/para  H2  converter2  at  T  =15  K  yields  a  flow  of  pre-cooled  pH2  gas  containing  ~  0.01  % 
residual  ortho-H2  and  ~  0.03  %  HD  (assuming  natural  isotopic  abundance)  and  only  sub-ppm 
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concentrations  of  condensable  impurities.  This  cold,  purified  pH,  gas  is  delivered  to  the  BaF2 
substrate  through  a  thin  wall  brass  tube  oriented  at  a  45°  angle  to  the  substrate  surface  normal, 
the  tube  terminates  about  3  cm  from  the  center  of  the  sample  substrate.  Because  of  the  sample- 
in-vacuum  design  of  our  cryostat,  during  a  deposition  the  cryogenic  thermal  isolation  vacuum 
must  be  maintained  below  ~  104  Torr  of  uncondensed  H,  gas  to  avoid  thermal  runaway.  An 
independent  flow  of  dopant  species  produced  by  any  of  a  variety  of  standard  methods  can  be 
readily  co-deposited  with  the  pH2  gas. 

The  majority  of  the  pH,  gas  freezes  to  form  a  curved  solid  deposit  on  the  "front"  side  of 
the  sample  substrate  and  holder,  depicted  schematically  in  Fig.  1 .  The  uncondensed  pH2  gas 
results  in  the  deposition  of  a  second  solid  H,  film  on  the  "back"  side  of  the  BaF2  substrate.  We 
will  show  below  that  for  our  aparatus  this  back  side  deposit  amounts  to  8.7  %  of  the  main  front 
side  sample  thickness,  and  that  this  11.5:1  front:back  thickness  ratio  is  roughly  independent  of 
the  pH,  inlet  flow  rate  over  the  1  to  60  mmol/hr  range.  The  amount  of  pH2  gas  which  “misses” 
the  sample  substrate  can  be  controlled  by  adjusting  the  position  of  the  brass  inlet  tube:  moving 
the  end  of  this  tube  closer  to  the  substrate  results  in  more  efficient  pH,  deposition,  but  also  in  a 
more  strongly  curved  front  side  deposit. 

The  IR  absorption  spectra  of  the  pH2  solids  used  to  calibrate  our  sample  thickness 
determinations  are  obtained  along  the  main  optical  axis  with  a  Fourier  transform  infrared  (FTIR) 
spectrometer  (Mattson  Nova  Cygni  120)  capable  of  0.09  cm  1  resolution,  which  simultaneously 
probes  both  front  and  back  side  pH,  samples.  We  employ  a  glowbar  source,  a  Ge  coated  KBr 
beamsplitter,  and  an  1N2  cooled  HgCdTe  detector  (16  pm  long-wavelength  cut-off)  to  obtain 
interfero grams  which  we  Fourier  transform  (triangle  apodization,  Mertz  phase  correction)  into 
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spectra  over  a  usable  range  of  800  to  5400  cm'1.  To  accommodate  the  IR  diagnostic  the  cryostat 
and  IR  optics  reside  inside  a  0.5  m3  polycarbonate  box  purged  with  a  constant  flow  of  dry  N2  gas. 

During  some  depositions  we  also  monitor  the  sample  growth  by  back-reflection 
interferometry  using  a  HeNe  laser  beam  incident  along  the  optical  axis  (9  =  0):43,44 

d  =  m  -  (8) 

2  n  cos(#) 

in  which:  d  is  the  pH2  sample  thickness,  m  is  the  number  of  fringes,  XHeNe  =  0.633  pm  is  the 
HeNe  laser  wavelength,  and  we  take  the  index  of  refraction  n  =  n^  =  1.14. 

We  also  present  recently  recorded  absorption  spectra  of  a  CO,/pH2  sample  obtained  with 
a  Bruker  IFS120HR  FTIR  spectrometer  using  a  Ge  coated  KBr  beamsplitter  and  two 
source/detector  combinations:  a  glowbar  source  and  a  different  1N2  cooled  HgCdTe  detector  (18 
pm  long-wavelength  cut-off),  and  a  tungsten  filament  source  combined  with  an  1N2  cooled  InSb 
detector. 


B.  Analysis  of  Spectrometric  Errors 

The  wedged  shape  of  the  front  side  solid  pH2  deposit  can,  in  principle,  introduce 
spectrometric  errors,  5A(v) ,  into  our  absorption  measurements:45 


8A(v)  «  -I-  [A(v)]2 

10  V  d  j 


(9) 


where  Ad  is  the  thickness  variation  across  the  sample  observation  volume;  however,  we  show 
here  that  under  our  experimental  conditions  such  errors  are  negligible.  We  assume  that  the  pH2 
flux  from  the  brass  inlet  tube,  and  hence  the  deposited  sample  thickness,  can  be  approximated  for 
small  angles  (j)  by:46 
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cos2  (<j>) 


(10) 


r" 

where  <j>  and  r  are  measured  from  the  centerline  and  the  open  end  of  the  inlet  tube,  respectively. 
The  front  side  pH2  deposit  depicted  in  Fig.  1  is  drawn  assuming  the  validity  of  Eqn.  (10)  and  a 
constant  pH2  sticking  probability;  the  shape  compares  well  with  visual  observations.  We  focus 
the  IR  probe  beam  to  a  3  mm  diameter  (or  smaller)  spot  at  the  sample,  yielding  an  upper  limit  of 
Ad/d  «  0.09  and  5A(v)  *  8x1 04  [A(v)  ]2,  which  we  disregard  for  A(v)  <  1. 

Our  use  of  photoconductive  HgCdTe  IR  detectors  raises  the  possibility  of  spectrometric 
errors  originating  from  (overdriven)  detector  nonlinearities.14'47’48  These  nonlinearities  primarily 
influence  the  strong  centerburst  region  of  the  interferogram,  which  Fourier  transforms  into 
distortions  of  the  spectral  baseline.  Indeed,  inspection  of  typical  single-beam  I0(v)  spectra  from 
both  the  Mattson  and  Bruker  FTIR  systems  reveals  false  non- zero  intensities  at  wavenumbers 
below  the  «  700  cm'1  BaF,  window  cut-off,  confirming  our  operation  of  these  detectors  in 
overdriven  IR  intensity  regimes.  This  problem  is  partially  mitigated  by  the  fact  that  in  any  given 
MIS  experiment  the  combination  of  pH2  sample  absorption  and  reduced  Fresnel  reflection  losses 
at  the  BaF2  substrate-to-vacuum  interfaces  modify  the  total  intensity  of  the  IR  beam  by  less  than 
±  10  %,  so  the  detector  must  follow  the  centerburst  signal  over  only  this  modest  range. 
Numerical  correction  for  the  resulting  baseline  distortions  is  complicated  by  the  partial 
compensation  for  detector  nonlinearities  which  occurs  during  the  phase  correction  stage  of  the 
FT  process.48  Instead  of  attempting  any  such  theoretical  corrections,  we  present  the  results  of 
additional  experiments  with  the  same  Mattson  FTIR  system  in  order  to  quantify  the  systematic 
errors  arising  from  detection  nonlinearity.  We  will  present  the  absorption  spectra  from  the 
Bruker  FTIR  system  without  performing  a  separate  independent  evaluation  of  nonlinear  errors. 
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To  this  end  we  measure  IR  absorption  spectra  of  a  nominal  5  %  CHC13/CC14  solution  held 
in  a  0  to  0.6  cm  variable  pathlength  cell  (Specac  7500)  equipped  with  non-rotating  KRS-5 
windows.  Single-beam  spectra  of  the  cell  filled  with  CC14  (obtained  at  each  pathlength) 
constitute  the  I0(v)  data.  We  calibrate  the  cell’s  length  scale  experimentally  from  the  pattern  of 
transmission  interference  fringes  in  IR  spectra  of  the  empty  cell  at  several  pathlengths,  d. 

d= _ m  - . —  (11) 

2ncos(0) [v,  -v2] 

in  which:  9  =  0  for  normal  incidence,  n  =  1 .0003  is  the  index  of  refraction  of  air,  and  m  is  the 
number  of  fringes  observed  over  the  vj  -  v2  wavenumber  range.  We  estimate  our  absolute 
pathlength  errors  as  ±  0.001  cm.  We  adjust  the  IR  intensity  at  the  HgCdTe  detector,  and  focus 
our  attention  to  the  4000  to  4500  cm 1  region,  in  order  to  reproduce  as  closely  as  possible  the 
same  overdriven  detector  conditions  as  encountered  during  the  pH2  experiments.  During  these 
experiments  the  CHC13  absorptions  and  mechanical  cell  replacement  errors  result  in  about  ±  10 
%  variations  in  the  centerburst  intensities  between  I0(v)  and  I  (v)  measurements. 

The  top  panel  of  Fig.  2  shows  a  correlation  plot  of  the  integrated  intensity  of  the  v,+v4 
CHC13  absorption49  vs.  sample  pathlength,  along  with  the  best  fit  linear  (dotted  trace)  and 
quadratic  (solid  trace)  polynomials: 

|A(v)dv:'  sy  =  a0  +  a1d  (12a) 

=  b0  +  b,d  +  b2d2  (12b) 

with  parameter  values:  a^  =  0.3142  cm’1,  a,  —  51.73  cm2,  b0  =  -0.0747  cm  ,  —  54.98  cm  ,  and 

b2  =  -5.041  cm'3.  We  note  that  for  d  =  0.573  cm,  the  peak  absorbance  is  A(4214  cm'1)  =  0.865,  so 
we  can  ultimately  generalize  our  results  from  pathlength  to  peak  absorbance  in  the  neighborhood 
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of  v  «  4200  cm'1  by  the  substitution:  d  =  1  cm  — »  A(v)  =1.51. 

At  first  glance  the  data  seem  to  fall  right  on  the  straight  line,  but  closer  inspection  reveals 
systematic  deviations  which  justify  the  use  of  a  quadratic  fit.  The  negative  value  of  the  b2 
coefficient  denotes  the  negative  curvature  consistent  with  a  mild  saturation  of  the  detection 
response  over  this  range  of  absorbances.  The  close  approach  by  the  fitted  quadratic  curve  to  the 
origin  and  the  uniform  distribution  of  the  residuals  in  the  bottom  panel  of  Fig.  2  confirm  that  this 
fit  adequately  describes  our  data. 

Recalling  Eqn.  (3),  the  quantity  of  interest  in  our  measurements  is  the  quotient  y/d,  which 
would  be  equal  to  the  slope,  y',  of  the  y  vs.  d  correlation  plot  in  the  absence  of  nonlinearities.  In 
practice,  the  slope  extracted  from  a  linear  fit  includes  an  error  which  we  can  estimate  using  the 
parameters  of  the  quadratic  fit.  The  slope  of  a  line  tangent  to  the  quadratic  fit  curve  at  any  given 
pathlength  is: 

y'  =  bj+2b2d  (13). 

By  numerical  comparison  we  find  that,  for  the  data  in  Fig.  2,  the  slope  of  the  best  straight  line  fit 
over  some  range  0  <  d  <  (fiax  is  approximately  equal  to  the  slope  of  the  line  tangent  to  the 
quadratic  fit  at  davg  =  d^/2.  Thus,  we  can  estimate  the  fractional  systematic  error  in  the  slope  of 
a  best  fit  straight  line  as  a  function  of  maximum  pathlength  or  peak  absorbance  as: 

-.0.09  <U  =  -0.06  A.„(v)  (14). 

y  b, 

As  expected,  the  fractional  error  in  the  fitted  slope  increases  for  thicker  samples  and  higher  peak 
absorbances.  We  will  use  Eqn.  (14)  below  to  estimate  the  magnitudes  of  systematic  nonlinear 
errors  in  our  pH2  experiments. 
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m.  RESULTS  AND  DISCUSSION 

A.  Single  and  Double  IR  Transitions 

In  pure  solid  pH,,  infrared  activity  in  the  homonuclear  diatomic  H,  molecules  arises  from 
dipole  moments  induced  by  intermolecular  interactions  with  neighboring  pH,  molecules.  An 
especially  helpful  and  succinct  summary  of  the  observed  IR  transitions  and  their  underlying 
physical  mechanisms  can  be  found  in  ref.  24.  Very  briefly:  the  single  and  double 
absorptions  are  zero-phonon  lines  involving  the  ro-vibrational  degrees  of  freedom  of  one  pH2 
molecule,  or  a  pair  of  pH,  molecules,  respectively.  Of  special  importance  here  is  the  fact  that  the 
interactions  responsible  for  single  transitions  are  subject  to  a  symmetry-dependent  cancellation 
effect;”  making  their  transition  intensities  very  sensitive  to  the  particular  structures  of  the  sites 
occupied  by  the  pH,  molecules.  In  contrast,  for  the  double  transitions  (in  which  a  pair  of  nearest 
neighbor  pH,  molecules  share  the  energy  from  a  single  photon)  the  transition  intensities  are  very 
insensitive  to  the  symmetry  of  the  surrounding  pH,  solid,  and  depend  primarily  on  the  solid  s 
bulk  density.  In  addition  to  these  zero-phonon  transitions,  broad  and  typically  strong  phonon 
branches  with  long  high-energy  tails  (“R-branches”)  are  also  observed. 

Fig.  3  shows  the  sharp  S,(0)  line  near  4486  cm'1,  and  the  associated  broad  SR(0)  phonon 
band  peaking  near  4550  cm’1,  as  well  as  the  two  double  transitions  which  are  the  main  focus  of 
this  manuscript:  the  Q,(0)+S0(0)  (4495  to  4520  cm'1)  and  S,(0)+S0(0)  (4825  to  4855  cm'1)  bands. 
These  double  absorption  features  are  not  spectrally  well-isolated  from  neighboring  bands,  so  we 
arbitrarily  choose  the  above  listed  wavenumber  ranges  as  convenient  band  integration  limits;  the 
corresponding  linear  integration  baselines  are  depicted  in  Fig.  3  as  dotted  line  segments.  We 
present  data  below  demonstrating  the  relative  merits  of  single  and  double  transitions  for  the 
determination  of  solid  pH,  sample  thicknesses. 
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Fig.  4  shows  the  effects  of  sample  annealing  on  the  intensity  of  the  S,(0)  feature  for 
another  rapid  vapor  deposited  pH2  solid.  During  the  temperature  cycle,  the  S,(0)  line  grows 
irreversibly,  roughly  doubling  in  integrated  intensity;  however,  the  integrated  intensity  of  the 
Q,(0)+S0(0)  band  (not  shown)  remains  constant  within  ±  2  %.  As  mentioned  above  in  the 
introduction,  we  have  previously  demonstrated1 '2'4,7'8  that  our  rapid  vapor  deposited  pH2  samples 
do  not  initially  solidify  into  the  thermodynamically  stable  hep  crystal  structure  obtained  by 
slowly  freezing  liquid  pH2.  Our  as-deposited  samples  are  neither  amorphous  nor  porous, 
however  any  close  packed  structure  such  as  hep,  or  fee,  or  conceivably  even  random-stacking  of 
close  packed  planes  (rep)  may  be  present.  We  have  also  previously  shown  that  annealing  of 
these  pH2  samples  to  above  T  «  4.5  K  irreversibly  converts  most  of  the  non-hep  regions  to  hep 
(however  —  1 00  ppm  concentrations  of  large  dopants  such  as  C60  can  inhibit  this  transformation3). 
The  presence  of  the  S^O)  line  in  Figs.  3  and  4(a)  requires  that  some  of  the  pH2  molecules  initially 
reside  in  sites  which  lack  a  center  of  inversion,21  Le.  in  hep  regions  or  perhaps  in  locally  “hep- 
like”  layers5’7  within  rep  regions.  The  irreversible  growth  of  the  S,(0)  line  shown  in  Fig.  4  is  thus 
consistent  with  the  permanent  conversion  of  some  or  all  of  the  non-hep  regions  into  hep  upon 
annealing.  We  observe  similar  annealing  behavior  for  the  U0(0)  line  near  1167  era'1.  Thus,  the 
variability  of  these  single  transitions  with  sample  microstructure  precludes  their  general  use  for 
determining  the  thicknesses  of  rapid  vapor  deposited  pH,  solids.  In  contrast  to  this  behavior,  the 
Q,(0)+S0(0)  and  S,(0)+S0(0)  bands  show  good  stability  over  the  life  of  a  ~  0.1  cm-thick  pH2 
sample,  seldom  varying  by  more  than  ±  5  %  even  after  annealing  at  T  =  4.8  K  for  1  hour 
(however  ~  10  urn-thick  samples  can  completely  sublime  away  under  these  conditions).  This  is 
due  to  the  fact  that  the  various  close  packed  structures  share  identical  nearest  neighbor 
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separations  and  bulk  densities. 

The  stability  of  the  intensities  of  the  double  transitions  is  further  illustrated  by  Fig.  5, 
which  shows  a  correlation  plot  of  the  integrated  intensities  of  the  Q,(O)+S0(O)  vs.  S1(O)+So(0) 
bands,  which  we  abbreviate  as  jAQiSdv  and  |AS+Sdv  ,  respectively.  We  note  that  Fig.  5  includes 
data  from  most  of  the  pH2  deposition  experiments  performed  in  our  laboratory  over  a  two-year 
period;  we  have  only  excluded  data  for  which:  JAQ+Sd  v  >  25  cm'1,  or  1AS+Sdi7  <  0.005  cm'1,  or 
from  spectra  obtained  at  resolutions  inferior  to  1  cm'1,  or  from  spectra  of  highly  scattering 
samples.  The  best  fit  straight  line  (constrained  to  pass  through  the  origin,  he.:  y  =  a,x)  has  a 
slope  of  a,  =  13.18  and  seems  to  represent  the  data  well,  even  for  the  thickest  samples. 

The  success  of  this  linear  fit  is  somewhat  surprising  as  we  might  reasonably  expect  the 
stronger  Q,(0)+S0(0)  band  to  exhibit  the  effects  of  nonlinear  detection  sensitivity  well  before  the 
weaker  S,(0)+S0(0)  band,  and  so  the  plot  might  reasonably  show  some  negative  curvature. 
Specifically,  the  thickest  samples  depicted  in  Fig.  5  exhibit  peak  absorbances  of  Amax(4510  cm'1) 
«  2.2  and  Amax(4843  cm'1)  *  0.17,  leading  to  expected  fractional  nonlinear  errors  from  Eqn.  (14) 
of  -  13  %  and  -  1  %,  respectively.  No  such  curvature  is  obvious  to  the  eye,  perhaps  being 
masked  by  the  scatter  in  the  data,  moreover  we  are  unable  to  detect  any  systematic  deviations 
from  linearity  via  numerical  analysis  of  these  data.  For  example,  restricting  the  data  sample  set 
further  to  JAQ+Sd  v  <  12  cm'1  and  then  to  |AQ+Sd  v  <  5  cm'1  yields  fitted  line  slopes  of  a,  =  12.86 
and  12.91,  respectively;  if  detection  nonlinearities  were  significant  we  would  expect  larger, 
rather  than  smaller,  a,  values. 

Fig.  6  shows  these  same  228  data  recast  as  a  histogram  of  the  quotient  jA^d  v  /[As+sd  v  . 
The  mean  value  is  p.  =  13.04  ±  0.08  (±  the  standard  error  of  the  mean50)  and  the  sample  standard 
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deviation  is  a  =  1 .22,  however  comparison  of  the  histogram  with  the  normal  distribution  having 
the  same  p  and  a  shows  that  the  measured  values  are  not  normally  distributed.  Most  of  the 
outliers  come  from  the  thinnest  samples;  for  example,  restricting  the  data  set  to  the  164  data  with 
p*  >  ^  cm'1  yields  p  —  13.02  db  0.06  and  o*  =  0.82,  and  eliminates  all  values  less  than  10  and 
greater  than  16.  Conversely,  restricting  the  data  set  to  the  64  data  with  jA^d  v  <  5  cm '  yields  p 
=  13.09  ±  0.24  and  a  =  1 .89,  so  even  the  data  from  these  thinnest  samples  show  roughly  the  same 
mean  value.  We  summarize  our  various  efforts  to  analyze  these  data  by  stating  that  our  best 
estimate  for  the  quotient  lAq^sd  v  /|As_sd  v  is  13.0  ±  0.4  (95%  confidence),  and  that  the  sample 
standard  deviation  is  <r  —  1 .2.  Since  our  analysis  of  these  data  do  not  reveal  any  effects  of 
detection  nonlinearity,  we  do  not  apply  any  correction  for  it. 

B.  Absolute  Sample  Thicknesses 

We  can  estimate  the  absolute  thicknesses  of  our  pH2  solids  from  the  pattern  of 
transmission  interference  fringes  appearing  in  some  IR  spectra  using  Eqn.  (11),  or  from  fringes 
recorded  by  back-reflection  HeNe  laser  interferometry  during  sample  deposition  using  Eqn.  (8), 
assuming  in  either  case  an  index  of  refraction  of  n  =  rp^  =  1.14.  Recalling  Fig.  1 ,  our  deposition 
method  forms  pH2  solids  on  both  the  “front”  and  “back”  sides  of  the  BaF,  substrate;  we  refer  to 
their  thicknesses  as  df-ront  and  dback,  respectively,  and  to  the  total  thickness  of  solid  pH,  in  the 
optical  path  as  =  dfront  +  dback. 

Fig.  7  shows  IR  spectra  for  three  different  pH2  solids;  the  intensities  of  the  pH2 
absorptions  in  the  4000  to  5000  cm'1  region  help  to  sort  out  the  evolution  of  the  pattern  of 
interference  fringes  with  increasing  total  sample  thickness.  For  very  thin  samples,  d^,  <  20  pm, 
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only  the  fringes  from  the  front  side  film  are  apparent,  as  shown  in  Fig.  7(a);  note  also  the  strong 
wavelength  dependent  scattering  losses  typical  of  slowly  deposited  pH2  solids.  Fig.  7(b)  shows  a 
somewhat  thicker  sample  which  exhibits  fringes  from  both  the  front  and  back  side  films;  for  such 
samples  with  20  jam  <  d,otaJ  <  70  pm,  we  can  simultaneously  determine  both  d^,  and  dback  directly 
from  the  IR  spectrum.  For  thicker  samples  the  curvature  of  the  front  side  deposit  increases  to  the 
point  that  its  transmission  interference  fringes  disappear,  leaving  just  the  fringes  from  the  back 
side  solid,  as  shown  in  Fig.  7(c).  Eventually,  for  even  thicker  samples  with  d,ota|  >0.1  cm,  the 
fringes  from  the  back  side  deposit  also  effectively  disappear,  leaving  a  conveniently  smooth 
spectral  baseline. 

Fig.  8  shows  the  back  reflection  HeNe  laser  interference  fringes  recorded  during  the 
deposition  of  a  Si/pH,  sample.  The  inset  shows  the  first  two  minutes  of  the  deposition;  during 
the  first  minute  fringes  due  to  the  growth  of  both  front  and  back  side  samples  are  apparent.  The 
direct  observation  of  fast  and  slow  fringes  in  the  back  reflection  measurements  confirms  our 
interpretation  that  the  interference  fringes  in  the  IR  spectra  arise  from  front  and  back  side  films. 

Fig.  9  shows  a  correlation  plot  of  the  integrated  absorbance  of  the  Q,(0)+S0(0)  transition 
vs.  the  interferometrically  determined  total  pH2  thickness,  for  all  1 1  samples  for  which  both  dfront 
and  dback  can  be  determined  simultaneously.  The  slope  of  the  best  fit  straight  line  through  the 
origin  is  a,  =  0.00821  cm' ‘/pm.  The  mean  of  the  quotient  jA<^sd  ?7cl0tal,  taken  over  these  1 1  data, 
is  p  =  0.00822  ±  0.00006  cm'Vpm;  the  sample  standard  deviation  is  cr  =  0.0002  cm'‘/pm.  For  the 
thickest  sample  depicted  in  Fig.  9,  the  peak  absorbance  is  ATnax(45 1 0  cm ')  ~  0.04,  resulting  in  a 
fractional  nonlinear  error  from  Eqn.  (14)  of  -  0.4  %,  or  -  0.00003  cm'Vpm.  Thus,  our 
recommended  value  is  lA^dF  /d,otal  =  0.0082  ±  0.0002  cm'Vpm,  or  equivalently,  82  ±  2  cm'2 
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(95%  confidence).  Finally,  incorporating  our  best  value  for  the  quotient  jAq+sd  v  /I A  s+sdv  we 
can  estimate  JAS+Sd  v  /d,ota|  =  6.3  ±  0.3  cm'2  (95%  confidence).  These  calculations  are 
summarized  in  Table  I. 

C.  Extrapolation  to  Thicker  pH2  Solids 

We  can  extrapolate  the  correlation  plot  of  Q,(0)+S0(0)  integrated  absorbance  vs.  pH2 
sample  thickness  up  to  nearly  0.1  cm-thick  samples  by  estimating  d,otal  from  our  interferometric 
measurements  of  dback.  We  obtain  the  ratio  of  the  ffont:back  sample  thicknesses  independently 
from  the  HeNe  back  reflection  measurements  and  from  the  IR  spectra  of  the  1 1  samples 
described  in  Fig.  9.  We  have  available  HeNe  back  reflection  interference  data  from  30  pH2 
depositions;  the  mean  of  4on/dback  over  these  data  is  u  =  1 1 .4  =  0.2,  with  a  sample  standard 
deviation  of  cr  =  0.8.  The  1 1  pH2  samples  for  which  the  IR  spectra  reveal  both  <4om  and  dback  yield 
a  mean  dfron/dback  of  li  =  12.3  ±  0.4,  with  a  sample  standard  deviation  of  a  =  1 .2.  By  trial  and 
error  we  find  that  a  value  of  dfron/dback  =  1 1.5  gives  values  of  d,otal  in  good  agreement  with  the 
extrapolation  of  the  |AQ+Sd  v  /dtQtal  =  82  cm"2  line,  as  shown  in  Fig.  10.  In  Fig.  10,  the  squares  are 
data  with  d,otal  estimated  as  d,otaJ  =  (12.5/1 1 .5)dfront,  and  the  triangles  show  values  of  estimated 
as  d,otal  =  12.5dback.  Of  course,  this  ffont:back  thickness  ratio  is  specific  to  this  particular 
experimental  setup.  More  important  is  the  apparent  absence  of  any  significant  curvature  for  solid 
pH2  thicknesses  approaching  0.1  cm. 

With  this  reassurance,  we  can  estimate  the  maximum  useful  thicknesses  over  which  the 
Q,(0)+S0(0)  and  S1(0)+S0(0)  double  transitions  can  be  used  to  determine  pH2  sample  thicknesses. 
If  we  adopt  the  criterion  A^  (v)  <  1  to  minimize  nonlinear  detection  errors,  then  we  can  use  the 
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Q,(0)+S0(0)  transition  for  d10ta]  up  to  about  0.13  cm-thick  pH2  solids,  and  the  S,(0)+S0(0) 
transition  for  up  to  about  1.7  cm-thick  solids: 

d,otal  (cm)  =  jA^gd  v  (cm'1)  /  82  (cm'2)  (15a) 

d,otal  (cm)  =  JAS+Sdv  (cm'1)  /  6.3  (cm'2)  (15b). 

Finally,  we  can  estimate  the  thickness  of  the  front  side  pH2  deposit  in  our  apparatus,  which  will 
contain  any  intentionally  introduced  dopant  species,  as: 

4.«  =  (11.5/12.5)  <U  (16). 


D.  Comparison  with  Literature  Values 

In  the  literature  on  the  spectroscopy  of  pure  molecular  solids,  two  alternative  forms  of  the 
integrated  transition  intensity  are  in  regular  use;25  retaining  symbols  consistent  with  our  previous 
definitions  we  have: 


^  2.303c  r  dv 

a  = -  A  (v)  — 

NpH2  d  Jband  v 


(17) 


in  which:  c  =  2.998x1 010  cm/s  is  the  speed  of  light,  and  a  has  units  of  cm3/s,  and: 

a=  f  a  (v)  dv  (18) 

£[  Jband 

in  which  a  has  units  of  cm'2.  This  integrated  transition  intensity  a  is  thus  related  to  the 
integrated  molecular  absorption  coefficient  defined  above  in  Eqn.  (3)  by: 


a=c*  La'(i7)  d"  =  y-  La>(7)  <19) 

m 

in  which  Vra  is  the  host  molar  volume  (Vra[pH2]  =  23.2  cm3/mol  at  lHe  temperatures,  ref.  18).  For 
reasonably  sharp  absorption  bands:25 
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a  « 


(20) 


ac 

NpH2  <  ?  > 

where  <  v  >  is  the  centroid  of  the  absorption  band.  For  the  Q,(0)+S0(0)  and  S,(0)+S0(0)  double 
transitions  of  interest  here  we  take  <  v  >Q+S  =  4507  cm'1  and  <  v  >s+s  =  4840  cm'1. 

Table  II  shows  a  summary  of  our  integrated  transition  intensity  calculations,  and  lists  the 
available  literature  experimental  and  theoretical  values  for  comparison.  Our  value  for  5q+s  and 
as+s  are  within  ±  25  %  of  the  previously  published  experimental  values,  however  these 
deviations  fall  well  outside  our  stated  95  %  confidence  limits.  We  note  that  the  experimental 
literature  values  were  obtained  from  spectra  recorded  using  pH,  sample  thickness  of  0.1  cm  for 
the  Q,(0)+S0(0)  transition,19  and  2.6  cm  for  the  S[(O)+S0(O)  transition.24  We  believe  that  our  use 
of  much  thinner  pH2  samples  helps  to  reduce  systematic  errors  due  to  detection  nonlinearities  and 
other  manifestations  of  the  finite  dynamic  range  of  spectrometric  apparatus. 

E.  Quantitative  Absorption  Spectroscopy 

We  conclude  our  discussion  by  presenting  step-by-step  sample  calculations  of  the  C  and 
0  atom  isotopic  abundances,  and  the  absolute  C02  dopant  concentration  derived  from  previously 
reported7  spectra  of  a  nominally  1 .2  ppm  C02/pH2  sample.  This  serves  as  a  definite  example  of 
the  methods  described  above,  documents  the  capabilities  of  our  experimental  apparatus,  and 
illustrates  certain  issues  peculiar  to  IR  absorption  spectroscopy  in  solid  pH2. 

Fig.  1 1  shows  the  v3  absorptions  of  the  12CI602  and  13C1602  isotopomers  at  0.008  cm'1 
FWHM  resolution,  not  shown  is  the  v3 12C160180  region  near  2327  cm'1  in  which  similar  features 
appear  with  a  peak-to-peak  S:N  ratio  of  about  3:1.  We  record  spectra  at  0.008  cm'1  FWHM 
resolution  using  a  glowbar  source  and  a  HgCdTe  detector  (40  co-added  scans  acquired  in  0.79 
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hr),  and  at  0.05  cm'1  FWHM  resolution  using  a  tungsten  filament  source  and  an  InSb  detector  (40 
co-added  scans  acquired  in  0.13  hr).  The  lower  resolution  data  obtained  with  the  more  sensitive 
InSb  detector  show  much  better  overall  S:N  ratios,  however  the  0.05  cm'1  resolution  is  not 
adequate  to  resolve  the  very  sharp7  C02  features.  Inadequate  spectral  resolution  can  be  an 
important  source  of  spectrometric  error  in  FTIR  measurements,14-47  and  is  an  especially  common 
pitfall  when  working  with  doped  pH,  solids. 

From  the  high  resolution  HgCdTe  spectra  we  extract:  lA(v312C1602)di7  =  0.1678  cm'1, 
jA(v313C1602)dv  =  0.0019  cm'1,  and  jA(v312Cl60180)dv  =  0.0005  cm'1;  using  3  cm'1 -wide 
integration  regions  of:  2343-2346  cm'1,  2326-2329  cm  *,  and  2277.5-2280.5  cm respectively. 
The  endpoints  for  linear  integration  baselines  are  obtained  by  averaging  over  1  cm'1 -wide  regions 
immediately  adjacent  to  the  integration  regions,  e.g.:  2342-2343  cm'1  and  2346-2347  cm 1  for  the 
12C1602  region.  We  estimate  the  95  %  confidence  limits  for  these  integrations  as  ±  0.0004  cm'1 
from  comparable  integrations  of  featureless  portions  of  the  2200  to  2400  cm'1  region.  These 
integrations  yield  spectroscopically  measured  isotopic  abundances  of:  12C  =  0.989  ±  0.004,  13C  = 
0.01 1  ±  0.002, 160  =  0.998  ±  0.003,  and  ,sO  =  0.003  ±  0.002,  in  excellent  agreement  with  the 
known51  natural  isotopic  abundances:  ,2C  =  0.9890(3),  ,3C  =  0.01 10(3),  160  =  0.9976(1),  and  180 
=  0.0020(1).  Thus,  the  natural  abundance  isotopomers  serve  as  convenient  internal  standards  for 
testing  the  importance  of  nonlinear  and  other  spectrometric  errors;  e.g.:  the  isotopic  abundances 
extracted  from  the  lower  resolution  InSb  data  show  much  larger  deviations  from  the  literature 
values. 

However,  the  lower  resolution  and  less  noisy  tungsten/InSb  data  are  preferred  for 
determining  the  thickness  of  the  pH2  solid  host.  Using  the  integration  limits  and  baselines 
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depicted  in  Fig.  3  we  extract  integrated  absorption  intensities  for  the  Q,(0)+S0(0)  and  S,(0)+S0(0) 
transitions  of:  lA^sdi7  =  11.79  cm'1  and  |As.sdv  =  0.932  cm'1,  respectively.  Applying  Eqn. 
(15a)  and  Eqn.  (15b)  yields  values  of  d.ota]  =  0.144  ±  0.004  cm  and  =  0.148  ±  0.007  cm, 
respectively.  As  these  thicknesses  are  a  little  larger  than  the  0.13  cm  maximum  we  recommend 
for  employing  the  Q,(0)+S„(0)  transition,  we  adopt  the  value  derived  from  the  S,(0)+S0(0) 
feature,  along  with  the  more  conservative  95  %  confidence  limits.  Finally,  applying  the 
apparatus-dependent  Eqn.  (16)  yields  the  thickness  of  the  doped  front-side  sample:  dfront  =  0.136 
±  0.007  cm. 

The  gas-phase  value52,53  of  the  integrated  molecular  absorption  coefficient  for  the  v3  mode 
of  C02  is  Jctv3 C02 ( v)  dv  =  5.49(0.21)xl07  cm/mol;  including  the  Polo-Wilson  factor  for  solid 

pEl2  from  Eqn.  (6)  yields  a  corrected  value  of  5.82(0.22)xl07  cm/mol.  Thus,  the  absolute  l2Ci602 
concentration  calculated  using  Eqn.  (3)  is  Cc02  =  4.9(0. 3)xl 0'8  mol/cm3  (95  %  confidence);  Eqn. 
(5)  yields  a  fractional  concentration  of  ,2CI602  =  1.13  ±  0.07  ppm,  in  good  agreement  with  the  1.2 
ppm  value  estimated  from  the  inlet  gas  quantities  of  0.1 1  pmol  C02  and  95  mmol  pH2.  The 
major  contributors  to  the  uncertainty  in  the  C02  concentration  are  the  sample  thickness  (5  %)  and 
the  literature  value  for  the  integrated  absorption  coefficient  (4  %). 

Similarly,  we  can  estimate  the  concentrations  of  the  other  C02  isotopomers  as:  I3C1602  = 
0.013  ±  0.003  ppm  and  12CI60180  =  0.003  ±  0.002  ppm.  Thus,  under  these  experimental 
conditions,  we  can  estimate  the  detection  limit  for  CO,  molecules  in  a  ~  0.1  cm-thick  pH2  solid 
as  ~  0.003  ppm,  with  S:N  «  3:1.  Employing  the  more  sensitive  InSb  detector  at  an  appropriately 
higher  resolution  would  extend  this  limit  to  below  ~  0.001  ppm,  which  is  comparable  to  the 
state-of-the-art  for  MIS  studies  on  doped  RGS  hosts  deposited  within  an  integrating  sphere.15 
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Combining  the  use  of  a  pH2  host  with  the  integrating  sphere  geometry  should  yield  even  lower 
detection  limits. 

IV.  SUMMARY 

Rapid  vapor  deposited  pH2  solids  offer  a  number  of  advantages  as  hosts  for  quantitative 
IR  absorption  spectroscopy.  In  particular,  the  sample  thickness  can  be  easily  determined  from 
the  observed  absorption  intensities  of  the  pH2  host  Q](0)+S0(0)  and  S,(0)+S0(0)  double 
transitions.  Tables  I  and  II  summarize  our  recommended  values  for  the  integrated  absorption 
intensities  of  these  two  bands,  which  we  calibrate  against  interfererometrically  determined 
sample  thicknesses.  Careful  evaluation  of  the  spectrometeric  errors  in  our  apparatus  encourages 
our  belief  that  our  new  values  are  more  reliable  than  those  published  previously.  Application  of 
Beer’s  Law  to  data  obtained  from  a  »  0. 1  cm-thick  C02/pH,  sample  demonstrates  our  ability  to 
detect  trapped  C02  molecules  at  concentrations  of  a  few  parts-per-billion  . 
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TABLES 


Table  I.  Summary  of  recommended  values.  We  abbreviate  Q,(0)+S0(0)  as  ;‘Q+S”  and 
Si(0)+S0(0)  as  “S+S.”  The  stated  uncertainties  indicate  the  95  %  confidence  limits. 


Ouantitv 

Recommended  Value 

Comments 

H2 

1.14 

refs.  38-40 

JAQ+Sd  v  //As^sd  v 

13.0  ±  0.4 

this  study 

jA^sdv/d,^ 

82  ±  2  cm'2 

this  study 

JAS+Sd  v  / d,otaJ 

6.3  ±  0.3  cm'2 

this  study 

^fron/^back 

11.5 

apparatus  specific 

Table  II.  Comparision  with  literature  values.  We  abbreviate  Ql(0)+S0(0)  as  “Q+S”  and 
S,(O)+S0(O)  as  “S+S.”  The  stated  uncertainties  indicate  the  95  %  confidence  limits. 


Ouantitv 

This  Study 

Literature 

Comments 

<  V  >Q+S 

4507  cm'1 

<  v  >S-b s 

4840  cm'1 

aS+S 

13.0  ±0.4 

^Q+S^  ^S+S 

14.0  ±0.4 

9.3 

expt.,  refs.  19, ' 

aQ+S 

1 89  ±  5  cm'2 

as+s 

14.5  ±  0.7  cm'2 

^Q+S 

4.84(±0.13)xl0'14  cm3/s 

4.2x1  O'14  cm3/s 

expt.,  ref.  19 

4.5x1 0'14  cm3/s 

theory,  ref.  20 

^S-fS 

0.35(±0.02)xl0'14  cmVs 

0.45x1 0'14  cm3/s 

expt.,  ref.  24 

0.97xl0'14  cm3/s 

theory,  ref.  24 

0.56x1 0'14  cm3/s 

theory,  ref.  26 
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FIGURE  CAPTIONS 


FIG.  1 .  Experimental  Diagram. 

FIG.  2.  Test  of  Mattson  FTIR  spectrometer  absorbance  linearity.  The  inset  shows  the  v,+v4 
CHC13  absorption  region  for  a  nominal  5  %  CHC13/CC14  solution  at  a  pathlength  of  0.573  cm. 
The  top  panel  is  a  correlation  plot  of  the  integrated  absorbance  of  this  feature  vs.  sample 
pathlength,  along  with  the  best  fit  (unweighted  least  squares)  straight  line  (dotted  trace)  and 
quadratic  polynomial  (solid  trace).  The  bottom  panel  shows  the  residuals  from  the  quadratic  fit. 
The  different  symbols  indicate  successive  experimental  runs. 

FIG.  3.  Condensed  phase  induced  features  in  the  IR  absorption  spectrum  of  a  vapor  deposited 
pH2  solid.  The  total  (front  +  back)  sample  thickness  is  dtotal  =  0.1 1  cm;  the  sample  was  produced 
in  a  deposition  time  of  tdep  =  0.5  hr  using  a  pH,  gas  flow  rate  of  ®pH2  =  125  mmol/hr.  The  entire 
spectrum  has  been  shifted  vertically  to  yield  zero  absorbance  at  4000  cm'1 .  The  dotted  line 
segments  indicate  the  baselines  used  for  integration  of  the  Q,(0)+S0(0)  (4495  to  4520  cm'1)  and 
S!(0)+S0(0)  (4825  to  4855  cm'1)  regions.  The  spectral  resolution  is  0.125  cm'1  FWHM. 

FIG.  4.  IR  absorption  spectra  of  a  d^  =  0.038  cm  thick  pH,  sample  showing  the  variability  of 
the  S,(0)  feature  with  sample  history.  Trace  (a)  is  for  the  as-deposited  sample  at  T  =  2.3  K,  trace 
(b)  is  the  same  sample  warmed  to  T  =  4.7  K,  trace  (c)  cooled  back  down  to  T  =  2.4  K.  The 
spectral  resolution  is  0.09  cm'1  FWHM. 
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FIG.  5.  Correlation  plots  of  the  integrated  intensities  of  the  Qi(0)+S0(0)  and  S1(0)+S0(0)  features. 
Shown  are  data  from  228  spectra  of  90  different  pure  and  doped  (Cdopam  <  0.1  %)  samples  at 
temperatures  between  2  and  5  K.  The  various  symbols  indicate  the  spectral  resolutions 
employed:  open  triangles  0.09  cm"1,  open  circles  0.125  cm'1,  closed  triangles  0.25  cm'1,  closed 
circles  1 .0  cm"'.  The  unweighted  least-squares  best  fit  straight  line  is  constrained  to  pass  through 
the  origin  (y  =  a,x)  and  has  a  slope  of  a,  =  13.18. 

FIG.  6.  Histogram  of  the  quotient  JAQ+Sd  v7jAs,sd  v  for  the  same  data  depicted  in  Fig.  5.  The 
solid  curve  is  the  normal  distribution  with  the  same  mean  (p  =  13.04)  and  standard  deviation  (cr 
=  1 .22)  as  the  data. 

FIG.  7.  IR  absorption  spectra  of  three  vapor  deposited  H,  solids  showing  the  evolution  of  the 
pattern  of  transmission  interference  fringes  with  increasing  sample  thickness.  Trace  (a)  is  for  an 
Al/pH2  sample  with  dfront  =  9.3  pm,  produced  in  tdep  =  1 .0  hr  using  a  pH2  gas  flow  rate  of  ®pH2  = 
0.9  mmol/hr.  Trace  (b)  is  for  a  B/pH,  sample  with  dfront  =  34.3  pm,  and  dback  =  3.0  pm,  tdep  =  1.0 
hr,  and  OpH2  =  2.6  mmol/hr.  Trace  (c)  is  for  a  Na/pH2  sample  with  dback  =  32.5  pm,  tdep  =  0.5  hr, 
and  OpH2  =  53  mmol/hr.  The  atomic  dopants  are  produced  by  pulsed  laser  ablation.  The  spectra 
are  presented  at  1  cm'1  FWHM  resolution. 

FIG.  8.  Back-reflection  HeNe  laser  interference  fringes  recorded  during  the  deposition  of  a 
Si/pH2  sample;  tdep  =  0.5  hr,  and  ®pK2  =  3 1  mmol/hr.  The  inset  shows  the  first  two  minutes  of  the 
deposition;  the  ratio  of  ffontback  fringes  is  roughly  11:1.  The  nearly  75  complete  back  side 
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fringes  indicate  that  dback  =  20.7  jam. 


FIG.  9.  Correlation  plot  of  lA^dy  vs.  dtotal  including  data  from  all  1 1  samples  for  which  both 
dfront  and  dback  can  be  determined  directly  from  the  transmission  interference  fringes  in  the  IR 
spectra.  The  unweighted  least-squares  best  fit  straight  line  is  constrained  to  pass  through  the 
origin  (y  =  a,x)  and  has  a  slope  of  a,  =  0.0082 1  cm'Vpm. 

FIG.  10.  Extrapolation  of  the  linear  fit  to  lAQ+sdv  vs.  d^.  The  circles  show  the  same  data  as  in 
Fig.  9,  for  which  both  d^,  and  dback  are  determined  simultaneously.  The  squares  show  values  of 
dt0taI  estimated  from  d^,  and  the  triangles  show  values  of  d,otal  estimated  from  dback,  assuming  a 
frontback  thickness  ratio  of  1 1 .5 : 1 . 

FIG.  11.  IR  absorption  spectra  of  a  natural  isotopic  abundance  C02/pH2  solid.  Panel  (a)  shows 
the  v3  transition  of  I2C1602  while  panel  (b)  shows  the  v3  transition  of  1jC1602.  The  sample 
consists  of  0.1 1  umol  C02  and  95  mmol  pH2  codeposited  over  a  0.5  hr  period.  The  spectral 
resolution  is  0.008  cm'1  FWHM. 
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Fig.  1.  Tam  and  Fajardo 
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